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Abstract. In this work, we calculate the contribution of dipole-dipole interactions to 
the optical nonlinearity of the two-dimensional random ensemble of nanoparticles that 
possess a set of exciton levels, for example, quantum dots. The analytical expressions 
for the contributions in the cases of TM and TE-polarized light waves propagating 
along the plane are obtained. It is shown that the optical nonlinearity, caused by the 
dipole-dipole interactions in the planar ensemble of the nanoparticles, is several times 
smaller than the similar nonlinearity of the bulk nanocomposite. This type of optical 
nonlinearity is expected to be observed at timescales much larger than the quantum 
dot exciton rise time. The proposed method may be applied to various types of the 
nanocomposite shapes. 
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In recent years, much attention has been given to optical nanostructures containing, 
for instance, semiconductor quantum dots. Frequently, such nanocomposites are 
organized as two-dimensional structures piMj. Normally, these structures are dense- 
packed systems of nanoparticles. The high concentration of nanoparticles will require 
taking into account the contribution of interparticle interactions, specifically long- 
range dipole-dipole ones, to the nanocomposite optical properties. Usually, the dipole- 
dipole coupling in the system of nanoparticles is treated by numerical methods |5l|6]. 
This approach is very popular for description of the systems of oscillating dipoles. 
Widely used approximations based on homogenization (effective medium theories) have 
restrictions on the size of inclusions [7], typical quantum dots do not fall in this range. 
In Ref. [8] mean-field theory was used in order to calculate the impact of dipole-dipole 
interactions in the materials doped with five-level nanoparticles. It is worth mentioning 
that the mean-field approaches are only suitable for the medium comprising a set of 
dipoles with some ordering and for the nonzero mean field. In this study, we obtain 
analytical expressions describing the effect of the dipolar interactions on the nonlinear 
dielectric susceptibility of quantum dot planar systems. 

In these nanocomposites, incident laser light can excite the electric dipole moment 
in the nanoparticles. Upon increasing the quantity of the dipolar excited nanoparticles, 
the dipole-dipole interactions in the nanocomposite become of particular importance. To 
calculate the contribution of the dipolar interactions to nonlinear optical characteristics 
of the two-dimensional system of randomly positioned cylindrical nanoparticles, the 
approach developed in Ref. [9j can be used. In this approximation, it is assumed that 
the mean value p of the induced electric dipole moment of the particle is proportional 
to the electric field amplitude Ein of the incident (external) radiation, 

P = EmOivEin, (l) 

where a is the dimensionless polarizability of the particle, v is its volume, Em is the 
dielectric function of ambient medium. Here p denotes the absolute value of the particle 
dipole moment p which can be arbitrarily oriented in three dimensions. This continuous 
approximation can be justified for the quantum dots possessing a number of exciton 
levels. The experimentally measured values of the polarizabilities of the CdSe quantum 
dots in terahertz or static electric fields lie within the range from 0.4 (Ref. [10]) to 3.6 
(Ref. [H]). The lifetime of the quantum dot excited state with dipole moment p is several 
orders greater than the period of an optical field oscillation. After the period much larger 
than the exciton rise time, typically ranging from tenths of picosecond [12l[l3] to tens 
of nanoseconds [H], a great number of the particles will be in the excited state and 
the system achieves a steady state. In the steady state, we suppose that the average 
quantity of the particles with the dipole moment does not vary and the system reaches 
the thermodynamic limit. Thus, on a timescale of the optical field oscillation, the 
polarized particles can be regarded as static dipoles. By this means, the quantum dots 
with the dipole moment will induce on the test particle, being located at the coordinate 
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where is the position vector of the l-th particle, is its dipole moment. In polar 
coordinates, has two components {ri,(j)i). Hereinafter, we use the orientation of the 
test particle dipole moment as the selected direction. As a first approximation, we 
assume that the local field E or induced nonlinearities do not change polarizability a of 
the individual particles. This field is random since the quantum dots are polarized 
arbitrarily and placed haphazardly. After obtaining the probability distribution 
function W{E) for the projection E of the field E onto the selected direction, the 
contribution of the dipole-dipole interactions to the nonlinear dielectric susceptibility 
of the nanocomposite can be calculated with the help of the statistical mechanics 
methods [9]. 

Let us consider the monolayer of cylindrical particles that lie randomly in a plane. 
These particles are assumed to have the electric dipole polarization ([1]) after exposition 
to laser light and after the transition to the steady state. Due to the probabilistic nature 
of the excitation, the number of the dipoles with one direction is balanced by dipoles with 
opposite alignment. Thus, the ensemble lacks the ordering and the mean of the random 
field E is zero. As shown in the Ref. [15] for such a system, W{E) becomes the Gaussian 
distribution at high nanoparticle surface concentrations c = Ns/S > 0.6, where N is 
the number of the particles having been excited, S is the area of the sample, s is the area 
of the particle. Hence, we cannot apply directly the approach of Ref. [9] exploiting the 
Gaussian distribution for the calculation of the third-order dielectric susceptibility as 
such high surface nanoparticle concentrations are hardly achievable in practice. In order 
to obtain W{E) at c < 0.6, making use of the negative cumulant expansion followed by 
the inverse Fourier transform was proposed |15] : 



where p is the Fourier-transform variable and the cumulants A„ are defined as follows 



r (p) is the distribution function for dipole moment orientations. This approach is useful 
for the case of the identical functions r (p) for all the particles. Integration over radius r 
in Eq. ([3]) begins from 2ro which is the minimal possible distance between two centers of 
the cylindrical particles having a radius tq. It should be emphasized that in the lack of 
the ordering in the system, cumulant expansion in Eq. (j2]) contains only even-numbered 
An [15] • Here, the absence of the ordering means that the net dipolar moment of the 
system is zero. 

The free energy density of the ensemble is found by the standard method: 
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where T is the temperature and ks is the Boltzmann constant, v = sh, h is the width 
of the monolayer (height of the particle). Projection P of the sample macroscopic 
polarization onto the selected direction for monochromatic radiation and non-absorbing 
medium can be obtained using the thermodynamic relation [16] 



P 



(5) 



where ()( denotes time averaging. It is worth recalling that in this model the dipole 
of the test particle can possess two opposite orientations with equal probabilities. To 
perform the time averaging, we check our results for the independence on the change of 
the sign of the exponential argument in Eq. H] (Ref. [9]). 

In the case of linear or circular light polarization, expanding P in a series 

P = XlEin + Xz\Ein?Ein + (6) 

it is possible to obtain cubic optical susceptibility of the system as well as higher 
order nonlinear susceptibilities. 

Upon differentiating ([5]) and expanding P in a series, we get 
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To simplify these expressions, let us invoke the Fourier integral representation of the 
Dirac 5-function [T7] 

5^"n^) = ^/(^P)"exp(zpP)dp (8) 
and the definition of the Dirac (5-function n-th derivative [17] 

/(P)(5(")(P)dE = (-l)"/(")(0). (9) 



Here / is some function. The use of the above two formulas reveals that integrals ([7]) 
obtained after the expansion have nonzero values if k = I and k = 0,2,4, .. ., namely 

Jq n = -27r, J2 ,2 = 4:71. 
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In the absence of the ordering, the first non- vanishing term in expansion ([S]) contains 
a third-order optical susceptibility X3(ci;; w, w, — w) (here u is the cychc frequency of 
the incident light). After simplifications we derive formula for the self-induced Kerr 
nonlinear susceptibility: 

It should be noticed that the same relationship for xa can be obtained using the Gaussian 
distribution as W{E). The next nonzero term in expansion (jH]) contains Xi'- 

X. = -^-^^^. (11) 

It should be noted that the above formulas are derived under condition Em > 0. In the 
above calculations we did not limit the number of the terms in the negative cumulant 
expansion in ([2]), so that the higher cumulants make contributions to higher-order 
susceptibilities. 

Further, let us consider some special dipole configurations that can occur if the 
monolayer of the quantum dots is integrated into a planar waveguide or a surface 
polariton propagates along an interface with the nanoparticles. When all the dipoles 
have in-plane alignment and they are coUinear with two equiprobable opposite directions, 
the distribution function r (p) for dipole moment orientations is given by 

r(7) = (12) 

where the angle 7 is measured from the selected direction which is determined by the 
dipole orientation. Then, 

A2 = 1 — • 13 

For the case of the dipoles lying at the interface between two half-spaces with ei and 

£2 [S], 

ei +£2 



2 ■ 

For the in-plane oriented dipoles having arbitrary directions. 



(14) 



^(7) = — , A2 = . (15) 

^ ' 277 ' 8(2ro)^5^ ^ ' 

Both of the above cases can occur when a TM-polarized wave propagates along the 
plane containing the nanoparticles. 

Another situation arises when TE-polarized wave propagates along the interface. 
In this case, the dipoles are expected to be oriented perpendicular to the plane. So we 
utilize function (fT2|) taking into account the fact that the polar axis is perpendicular to 
the plane. Then, 
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For the medium consisting of two half-spaces with ei and 62 jlSj, 



As one can see, all of these variants of A2 differ little. 

By way of illustration, in the case of the dipoles with in-plain collinear polarizations, 
substituting v = sh, s = vrrg, and ( IT3l) in ( ITOl) we get 



It is obvious that the obtained expression for the Kerr nonlinear susceptibility does not 
depend explicitly on the particle radius. However, one should keep in mind that the 
polarizability of the particles strongly depends on their size. Provided c = 0.05, a = 0.4, 
h = 5 nm, Em = 2.3, T = 300 K using Eq. (fT8|) we arrive at ^3 ^ — 8 x 10^^^ esu. This 
value is of the same order as observed by experiment with nanosecond laser pulses for 
cadmium chalcogenide nanocomposites (— 10~^^ — 10~^° esu) [T9ll20] . 

As well as for the bulk nanocomposite, the contribution of the dipole-dipole 
interactions to the Kerr optical nonlinear susceptibility of the planar system has a 
negative sign. This is due to the disarranging effect of the random field of the 
induced dipoles. It is worth mentioning that the dipolar optical nonlinearity would 
become significant at the time intervals which exceed the exciton rise time many- 
fold. Comparison with the results of Ref. |9] shows that for the same r and h the 
contribution of the dipole-dipole interactions to the Kerr susceptibility of the two- 
dimensional nanocomposite is several times less than for the bulk nanocomposites. 
This is especially important for the nanoparticles which acquire the large values of 
the dipole moment, for instance, the II- VI semiconductor quantum dots [2^. Often, 
under femtosecond laser pulses, these quantum dots exhibit the positive real part of the 
third-order susceptibility [22l|23], but at nanosecond pulse durations they show negative 
nonlinear refractive index [T9l[20]. This phenomenon may be attributed to the effect of 
the dipole-dipole interactions. 

If the dipolar optical nonlinearity is undesirable, the organization of the 
nanocomposites, which contain the nanoparticles with the great magnitudes of the 
induced dipole moment, as planar structures is justified. However, one would take 
into account that other types of optical nonlinearity also may decrease with reducing 
the number of dimensions. 

As one can see from Eq. ffTOj) . X3 is proportional to c^. It is the dependence on 
concentration that is the key feature of the contribution of the interparticle dipole- 
dipole interactions to the optical nonlinearity of the composite. This is reasonable since 
the number of interparticle couplings varies as the square of the quantity of the particles 
in the sample. As one could expect, other types of optical nonlinearities would linearly 
depend on the nanoparticle concentration. 

It should be underlined that the proposed here approach to the calculation of the 
dipole-dipole interaction contribution to the nonlinear susceptibility of nanocomposite 
can be applied not only to the two-dimensional geometry but also to other configurations. 



X3 = 



USkeT 
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The method, based on the cumulant expansion of the characteristic function, can be 
also utihzed for other types of interparticle interactions, such as the quadrupolar ones, 
since there is no need to determine probabihty distribution function W{E) exphcitly. 

In conclusion, the analytical expressions for Kerr optical nonlinearity ^3 of the 
two-dimensional system of the randomly located quantum dots have been derived. 
It has been shown that for the different types of the light polarization, Xd. changes 
insignificantly. The dipole-induced optical nonlinearity of the two-dimensional sample 
is several-fold smaller than for the bulk one so that this can be used in practice. 
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